function bnet = mk_cancer()
N=5;
dag = zeros(5,5);
dag(1,2)=1;
dag(1,3)=1;
dag(3,4)=1;
dag(2,4)=1;
dag(3,5)=1;
node_sizes= zeros(1,5);
node_sizes(1)=2;
node_sizes(2)=2;
node_sizes(3)=2;
node_sizes(4)=2;
node_sizes(5)=2;
bnet = mk_bnet(dag, node_sizes);
CPT= zeros(2,1);
CPT(1)=0.2;
CPT(2)=0.8;
CPT=reshape(CPT, 1, prod(size(CPT)));
bnet.CPD{1} = tabular_CPD(bnet, 1, 'CPT', CPT);
CPT= zeros(2, 2);
CPT(1,1)=0.8;
CPT(2,2)=0.8;
CPT(1,2)=0.2;
CPT(2,1)=0.2;
CPT=reshape(CPT, 1, prod(size(CPT)));
bnet.CPD{2} = tabular_CPD(bnet, 2, 'CPT', CPT);
CPT= zeros(2, 2);
CPT(1,2)=0.8;
CPT(1,1)=0.2;
CPT(2,2)=0.95;
CPT(2,1)=0.05;
CPT=reshape(CPT, 1, prod(size(CPT)));
bnet.CPD{3} = tabular_CPD(bnet, 3, 'CPT', CPT);
CPT= zeros(2, 2, 2);
CPT(1,2,1)=0.8;
CPT(1,2,2)=0.2;
CPT(2,1,2)=0.2;
CPT(2,2,1)=0.05;
CPT(1,1,1)=0.8;
CPT(2,1,1)=0.8;
CPT(1,1,2)=0.2;
CPT(2,2,2)=0.95;
CPT=reshape(CPT, 1, prod(size(CPT)));
bnet.CPD{4} = tabular_CPD(bnet, 4, 'CPT', CPT);
CPT= zeros(2, 2);
CPT(1,1)=0.8;
CPT(2,2)=0.4;
CPT(1,2)=0.2;
CPT(2,1)=0.6;
CPT=reshape(CPT, 1, prod(size(CPT)));
bnet.CPD{5} = tabular_CPD(bnet, 5, 'CPT', CPT);
